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Abstract 



There is a good deal of current interest in the condensate < A^Aa > which has been seen 



. to play an important role in calculations which make use of the operator product expansion. 

That development has led to the publication of a large number of papers which discuss how that 

P-i! 

condensate could play a role in a gauge- invariant formulation. In the present work we consider 

<D ' 

■ gluon propagation in the presence of such a condensate which we assume to be present in the 

vacuum. We show that the gluon propagator has no on-mass-shell pole and, therefore, a gluon 

X 

cannot propagate over extended distances. That is, the gluon is a nonpropagating mode in the 



gluon condensate. In the present work we discuss the properties of both the Euclidean-space and 
Minkowski-space gluon propagator. In the case of the Euclidean-space propagator we can make 
contact with the results of QCD lattice calculations of the propagator in the Landau gauge. With 
an appropriate choice of normalization constants, we present a unified representation of the gluon 
propagator that describes both the Minkowski-space and Euclidean-space dynamics in which the 
< A^Aa > condensate plays an important role. 

PACS numbers: 12.38.Aw, 12.38.Lg 
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I. INTRODUCTION 

Recently studies making use of the operator product expansion (OPE) have provided 
evidence for the importance of the condensate < A a ^A% > [1-3]. (There is a suggestion that 
such a condensate may be related to the presence of instantons in the vacuum [4].) The 
importance of that condensate raises the question of gauge invariance and there are now a 
large number of papers that address that and related issues [5-19]. We will not attempt to 
review that large body of literature, but will consider how the presence of an < A a ^A^ > 
condensate modifies the gluon propagator and the vacuum polarization function in QCD. We 
may mention the work of Kondo [7] who was responsible for introducing a BRST-invariant 
condensate of dimension two, 



where c(x) and c(x) are Faddeev- Popov ghosts, a is the gauge-fixing parameter and fl is the 
integration volume. Kondo points out that fl reduces to A 2 min in the Landau gauge, a — 0. 
Here the minimum value of the integrated squared potential is A 2 miw which has a definite 
physical meaning [7]. 

For recent discussion of the role of various vacuum condensates in QCD one may refer 
to Refs. [20] and [21]. (In these works the value given for the gluon condensate is 
< (a s /ir)G 2 >= 0.009 ± 0.007 GeV 4 .) In our early work [22] we assumed that the gluon 
condensate carried little or zero momentum. The vector potential of the theory was divided 
into a condensate field, A"(x), and a fluctuating field, A®(x). The field A.^(x) is independent 
of x and has zero vacuum expectation value in our model.. 

We define an order parameter, in a covariant gauge: 




(1.1) 



< vac | A"(0)A*(0) I vac > 



8 4 




(1.2) 



The field tensor for QCD is given by 



G%(x) = d,A a v {x) - dvAltx) + gf abc Al(x)A:(x). 



(1.3) 



We insert 



A«(s) = A« + ^(a;) 



(1.4) 
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into Eq. (1.3) and define 



G%,(x)=G% + g^(x) 1 (1.5) 

where 

G% = gf ab XA c u (1.6) 

is the condensate field tensor. We stress that, if the zero- momentum mode is macroscopically 
occupied, A^ and G^ u may be treated as classical fields. However, we must maintain global 
color symmetry and Lorentz invariance when using such fields. 

As noted above, in our model [22], the gauge-invariant condensate parameter 

< vac |: g 2 {jl 2 )G a ^ v {Q)G^{Q) :| vac >, is related to the condensate parameter, 

< vac |: g 2 (/i 2 )A^(0)A^(0) :| vac >. This relation follows from our assumption that the 
condensate is in a zero-momentum mode. Since we have a phenomenological value for 

< vac |: g 2 {jx 2 )G a {ti)G% v {<S) :| vac >, obtained from QCD sum-rule studies [23], we can 
obtain a value for < vac |: g 2 (/2 2 ) A^O) A%(0) :| vac > by the following procedure. 

Using the assumption that the condensate carries zero momentum, we identify the con- 
densate contribution as 

' < vac\:g 2 (f, 2 )G a (0)G^(0):\vac> 



4?r 2 ^' ^ 

= JL < vac |: ^ 2 (/i 2 )G^0)G^0) :| vac > (1.7) 

= 4^2 J2f abC f ab ' C ' < vac I SVK(0K(0)A£(0)A C ^(0) | V ac > . (1.8) 
We may write 

< vac | AJ(0)A£(0)A£(0)AJ(0) | vac > 

= j^^^hA'c + g/g v v h b i5cc> + g» v g/5 bc ,6 cb ,}. (1.9) 

We have previously calculated matrix elements of this type by several methods. In one 
work we calculated matrix elements of the condensate potential after constructing | vac > 
as a coherent-state in the temporal gauge [24]. In another work [22] we wrote A^(0) = 
0o^j where r)°r)% — — 1. In the latter scheme 77® was averaged over the gauge group when 
calculating matrix elements of products of condensate fields. (One way to check the factor 
(32)(34), which appears in the denominator of Eq. (1.9), is to set b = c, ji = u, p = rj 



and b' = c' and sum over identical indices.) We may insert the vacuum state between the 
operators to obtain 

< vac | Aj(0)A£(0) | vac >< vac | A£,(0)Aj'(0) | vac >= <f%. (1.10) 

This then agrees with the result obtained when evaluating the right-hand side of Eq. (1.9). 
Now, using Eq. (1.9) in Eq. (1.8), we find 

^ < vac |: sV)G^(0)Gf (0) :| vac >= (4)( ^ (yV)0§) 2 , U-U) 
from which we obtain 

<7 2 (/2 2 )tfb = 1.34(GeV) 2 . (1.12) 

(Here we use the renormalization point jj 2 ~ 1 GeV 2 .) We will make use of these results in 
the following. 

In this work we discuss the form of the gluon propagator in some detail. We also contrast 
the structure of the propagator in QCD and QED. In this comparison the distinction between 
theories with and without boson condensates is particularly clear. A characteristic of a 
theory with condensates is the appearance of a term proportional to g^ v 5 A {k) in the gluon 
propagator. This term describes the macroscopic occupation of the zero-momentum mode 
and provides a covariant representation of the effect of the condensate in modifying the 
structure of the propagator. 

The organization of our work is as follows. In Section II we review the introduction of 
the vacuum polarization tensor in the case of QED. In Section III we discuss the vacuum 
polarization tensor for QCD and in Section IV we review the Schwinger mechanism for 
dynamical mass generation for gauge fields [25]. In Section V we define a dielectric function 
for QCD and present the results of our calculation of that quantity. In Section VI we 
provide values of the gluon propagator in both Euclidean and Minkowski space and make 
some comparison to the propagator obtained in lattice simulations of QCD. Finally, Section 
VII contains some further discussion and conclusions. 
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II. THE PHOTON PROPAGATOR AND THE DIELECTRIC FUNCTION IN 
QED 



In this section, we review standard results for the photon propagator in QED. (This 
material is available in the standard textbooks.) The propagator may be written as 



iV™{k) = -i 



(fiTtf k^k u /k 2 ) \k^k v 



k*(l-ni m (k*)) (k? + ie)\ 
where Ul m (k 2 ) is a finite quantity with the following limits (to order a), 

a k 2 



n em0 2 ) 



157T m? 



k 



0; 



a 



In 



3ir \m 

It is useful to define the dielectric function, 



5a 
9tt' 



k — > — oo. 



so that 



n{k 2 



<k 2 ) = [1 - TlUk 2 )} 



a k 2 

157T m z 



0, 



and 



K{k 2 ) 



l-^ln / " 



ha 

+ 57 + 



3/T \ m z 

A charge placed in the vacuum gives rise to a potential, 

e 



k 



— oo. 



V(k) = 



K(k 2 ) | k 



which, for small distances, behaves as 

V(k) 



k I 2 



(2.1) 



(2.2) 
(2.3) 



(2.4) 



(2.5) 



(2.6) 



(2.7) 



(2.8) 



This is the standard result, which indicates that QED becomes strongly coupled at short 
distances. We make the observation that D e ^{k) has a pole at k 2 =0. We make this ap- 
parently trivial observation, since we wish to demonstrate that in our model there is no 
corresponding pole in the gluon propagator in QCD - that is, the gluon becomes massive 
via the Schwinger mechanism [25]. 



For completeness, we note that the polarization tensor has the form, 



A' 2 



and we have 



<DP(*) = iVl{k) + [<n^(fc)] iV£(k), 



1)V 



where 



iV°(k) = -i 



k 2 + it 

It is also useful to rewrite Eq. (2.10) as 

W ltv (k)=iD^(k)+iV^(k) 



i^Qixv k^k v /k 2 ) \k^k v 



(k 2 + ief 



iU(k 



iV$(k). 



l-U(k)/k 2 _ 

The term [—i\k^k u /(k 2 + it) 2 } is common to both sides. Therefore we may put 

vlM = -[g»v - k»k v /k 2 \v T {k), 



(2.9) 



(2.10) 



(2.11) 



(2.12) 



(2.13) 



V^{k) = -[g, v -k,K/k 2 \V^\k), 



with 



V T (k) = 



k 2 -U(k)+ie' 



and 



Vf{k) 



k 2 + ie 



to obtain the relation 

U(k 2 ) 



[V?(k)}^[V T (k)-V^(k)][V^(k)]'\ 



(0), 



(2.14) 



(2.15) 



2-16) 



(2.17) 



l-U(k 2 )/k 

The left-hand side of Eq. (2.17) is related to the time-ordered product of the currents, 



1^ 



U(k 2 



l-U(k 2 )/k 2 



i I d 4 xe ik ' x < vac | T[f(x)f '(0)] | vac> . (2.18) 



Note that n(fc 2 ) is the irreducible self-energy, while the matrix element of the time-ordered 
product of the currents gives rise to a reducible form. 
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We remark that the equation for the vector potential is 

A v {x) = -f{x). (2.19) 



-g^D + (1 - \)d»d u 
A 



Thus 



d,f{x) = ±U(d v A v {x)). (2.20) 

However, in QED we have 

d^fix) = (2.21) 
as an operator relation from which it follows that 

□(^"(a;)) = (2.22) 

is an operator relation. 

Since the current is explicitly conserved in QED, Tl(k 2 ) is independent of the gauge-fixing 
parameter. In general, we can write 

f{x)=f T {x)+f L {x), (2.23) 

where 

j£( x ) = (g^-d^dnUx), (2.24) 

and 

f L (x) = d^d^Ux). (2.25) 

We see from Eq. (2.21), that f L (x) = in QED. 

In the case of QCD the situation is more complicated since current conservation does not 
appear at the operator level in a covariant gauge. Gauge fixing breaks the general local gauge 
invariance of the theory. (However, gauge fixing is necessary for covariant quantization, since 
without such a procedure one finds that the momentum conjugate to Aq(x) is zero.) In 
QCD one usually uses path-integral quantization. That formalism leads to the introduction 
of ghost fields. These fields insure unitarity for the gluon channels. 

Here we will follow Lavelle and Schaden [26] and define the nonperturbative gluon prop- 
agator 

Vnonpert^ = _ ^(fc 2 ). (2.26) 
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The nonperturbative propagator is transverse in any covariant gauge. That is a consequence 
of the Slavnov- Taylor identities which state that the full and the perturbative longitudinal 
propagators are the same in any covariant gauge. Therefore, the difference appearing in 
Eq. (2.26) yields a purely transverse result. In the following discussion we will consider 
the calculation of the nonperturbative gluon propagator and insure that our propagator is 
transverse. 

III. COVARIANT QUANTIZATION IN QCD 

Consider the Yang-Mills Lagrangian for a SU (3) color theory without quarks. We have 

c(x) = -\g; u (x)g^(x) - ^[d,A»{x)f - dM^Mx) 

+gf a bc[dnMx)Mx)A»(x)], (3.1) 

where (j> a (x) and 4> a (x) are ghost fields and A is a gauge- fixing parameter. Now with 

Gf(x) = PA&x) - d»A»{x) + gf abc A^x)A:(x), (3.2) 

and using Eq. (3.1), we have 

d,G^{x) + \d»{d»A»{x)) = f(x) (3.3) 

= -gf^{x)G^{x) -gr bc [d»4> b (x)}cf> c (x), 

which we write as 

-<rn + (i - A a u (x) = -r a {x), (3.4) 

where 

J»{x) = -gf abc A b v {x)G^{x) - gf abc d u (A^x)A^x)) - gf^M^M*)- (3-5) 
Note that J£(x) is conserved in the classical theory. We write 

J£(*) = JW + J W ( 3 - 6 ) 
using definitions analogous to those in Eqs. (2.24) and (2.25). Then 



d.J^x) = V£» 
8 



(3.7) 



and 



<VM = ja(d^(x)). 



(3.8) 



As in QED, the constraints imposed by current conservation in the physical Hilbert space 
can be maintained by calculating with (x) instead of J^(x). (We remark that the ghost 



fields do not contribute to Jj<^{ x )-) Note that 



(3.9) 



A 



nduAUx)) = r L (x), 



and 



-U&AlJx)) = d^x), 



which also follows from Eqs. (3.7) and (3.8). 
We may also write 

U(k 2 ) 



[g^ - k^jk 2 " 



l-U{k 2 )/k 2 



S ab = [g^k 2 - k»k 



Mk 2 ) 



ab 



l-U^k 2 ) 

= i J d 4 xe ihx < vac | T[J» a (x)J» 

where the polarization tensor is defined as 

Kb(k) = (g^-k^/k 2 )U(k 2 )5 ab , 

= (g^k 2 -k^)U l {k 2 )5ab. 



(3.10) 



(3.11) 



(3.12) 

vac >, 

(3.13) 
(3.14) 



The division of the field into transverse and longitudinal parts does not have the same 
utility in QCD as in QED, since JT t a( x ) anc ^ ^L,a( x ) have a nonlinear dependence on the 
gluon field. Therefore, the field equations do not separate into transverse and longitudinal 
equations in the case of QCD. The quantities Tl(k 2 ) and IIi (A; 2 ) are defined in terms of 
conserved currents, J^ a (x). If we work with Ja{x) rather than Jj, a (x), the ghost fields will 
insure that U.^(k) has a transverse structure. However, if one does not insure the constraint 
< tp | d^A^x) | if; > — 0, one finds a dependence on the parameter A in U(k 2 ). While 
the presence of ghosts insure unitarity relations, they do not serve to impose the constraint 
d^A^{x) = 0. Our calculation corresponds to a diagrammatic analysis, made in the Landau 
gauge, with condensate ghosts added to insure the transverse nature of U!^(k) [26]. 
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IV. DYNAMIC MASS GENERATION VIA THE SCHWINGER MECHANISM 



The term in J„(x), — g 2 f abc f cde A b (x)A^(x)A^ (x), will give rise to a gluon mass term 
if there is a gluon condensate in the QCD ground state. We find a contribution to the 
conserved (transverse) current of the form 



J£(x) = -ml 



y □ 



•A u (x), 



where 



m 



G 



32 



(4.1) 



(4.2) 



This corresponds to a contribution to the polarization tensor of the form 



nl(*0 = 6 ab (g^ - k^k u /k 2 )m 2 G , 



(4.3) 



with uiq = 614 MeV, if we use Eq. (1.12). [See Fig. 1]. It is, therefore, useful to define 



n A (fc 2 ) 

k 2 



where the second term does not have a pole as k 2 — > 0. We also have 



U(k 2 ) =m 2 G + U A (k 2 ). 



(4.4) 



(4.5) 



The appearance of a pole at k 2 = in TIi (fc 2 ) defines the Schwinger mechanism [25]. 
It is also useful to subtract the quantity given in Eq. (4.1) from the current and define 



(4.6) 



in momentum space. Then we have 



{k 2 - m 2 G ){g^ - Wk u /k 2 ) + -k"k v 

A 



K{k) = -#(*). 



We now write a first-order propagator as 



k 2 — m 2 G 



+ A 



WW 



(k 2 + ieY 



Sab, 



and also write 



iV{k 2 ) = iV {1 \k 2 ) + iV {1 \k 2 ) 



in A {k 2 ) 
L _ n A (fe 2 ) 

k 2 —m z G 



iVW(k 2 ), 



(4.7) 



(4.8) 



(4.9) 



10 




FIG. 1: a) Calculation of the gluon self-energy in the condensate. The first term shows the origin 
of the gluon mass term in the mean-field approximation. The dashed line refers to a condensate 
gluon of zero momentum. In the second part of a) we show a contribution to the (irreducible) 
polarization tensor in the single (condensate) loop approximation, b) Diagrams which contribute 
to the polarization tensor in QCD. The wavy line is a gluon, the solid line is a quark, and the third 
diagram represents the gluon field, c) Some corrections to the diagrams of (b) due to the presence 
of a gluon condensate. 
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where 

vW(k 2 ) = -[k 2 -m 2 G ]-\ (4.10) 

and 

V(k 2 ) = -[k 2 -m 2 G -U A (k 2 )]- 1 . (4.11) 

Thus we see that, after absorbing the mass term in the quantity TlA(k 2 ) is related 
to the time-ordered product of the Jt^( x )'- 

5 ab = i J d 4 xe ik - x < vac | T[J^ a (x)J^ b (0)} | vac > .(4.12) 

Equation (4.12) is a generalization of Eq. (2.18) and reflects the presence of a condensate in 
the QCD vacuum which makes the gluon massive. We also remark that, after gauge fixing, 
the theory with ghost fields has a form of gauge invariance - the BRST gauge symmetry. 
This symmetry allows one to derive the analog of the QCD Ward identities in QCD - the 
Slavnov- Taylor identities. (We note that the ghost condensate introduced in Ref. [26] is 
BRST invariant.) 

V. THE GLUON PROPAGATOR AND THE DIELECTRIC FUNCTION IN QCD 

One is tempted to write the analog of Eq. (2.1) in the case of QCD. However, if there 
is a gluon condensate present, there is an essential modification to be considered. We recall 
that we found it useful to divide A a {x) into a condensate field, A" (re), and a fluctuating 
field, A a „{x). (We made the assumption that the condensate field is in the zero-momentum 
mode and therefore, A£(x) is independent of x.) 

Thus, in coordinate space, we have 

iVf u (x,x') = <vac\ T[A«(x)A b v (rf)\ \ vac > (5.1) 
= < vac | A°A* \vac> + < vac | T[A a ^(x)A b u (x')} | vac > (5.2) 
= - 9 -f^f+ < vac | T[Al{x)A h v {x')\ | vac > . (5.3) 

Our expression for the gluon propagator in momentum space is then 

^V;l(k) = - 9 -f ( f ) 2 6 f(2 7 rr5^(k)-^ 



[g^ - k»k v /k 2 



Tl A (k 2 



1 - 



n A (fc 2 ) 



\Qim> k^k u / k ) 
A; 2 [1 -IIi (A; 2 )] ' '\k 2 + ie) 



+ A- 



kp}z v 



(5.4) 
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We see the first characteristic difference when we compare Eq. (5.4) with Eq. (2.1), that 
is, the presence of a delta function. Whether such a term is present depends on whether or 
not one has a ground-state condensate in the zero-momentum mode. 
We can define a QCD dielectric function: 

k(A; 2 ) = [1-IIi(A; 2 )]. (5.5) 

As noted earlier, the Schwinger mechanism refers to the a fact that, if n(k 2 ) has a pole at 
k 2 = 0, the gluon has a dynamical mass and the pole at k 2 = in V^ v {k) disappears. In an 
earlier work we found that m 2 G = (9/32)g 2 0Q. As we will see in this work 



:>. , ">r, W 



K(k 2 ) = 



k 2 k 2 — m 2 G 



(5.6) 



where r] 2 = (3/32)^ 2 0g. (Thus m 2 G = 3r] 2 .) The result given in Eq. (5.6) follows from the 
calculation of the diagrams in Fig. 1 subject to the constraints required in the covariant 
formalism. The quantity n(k 2 ) defined in Eq. (5.6) is shown in Fig. 2. It is interesting 
to note that if the Schwinger mechanism is operative, that is, if there is a pole in K,(k 2 ) at 
k 2 = 0, one needs an additional singularity to avoid having a zero in nik 2 ). Such a zero would 
imply that gluons could go on-mass-shell, a clearly unsatisfactory result. It is gratifying that 
at the next level of approximation (one condensate loop) one finds the necessary singular 
term that maintains the relation n(k 2 ) ^ 0. 

We do not have all the terms contributing to K(k 2 ). For example, there will be terms of 
order g 2 {(f^/k 2 ) or g 2 \n{(f)^/ \ k 2 |) in the deep Euclidean region. The origin of such terms 
may be seen in Fig. 1, where we have shown how the presence of the condensate can lead 
to (power) corrections to the asymptotic behavior of the polarization tensor in the region 
k 2 — ► — oo. 

We also note that the only way to form a small parameter in this model is to construct 
the ratio [g 2 4>1/{— k 2 )} which is small for large spacelike k 2 . The nonperturbative analysis 
is clearly not an expansion in a small parameter. That is characteristic of nonperturbative 
approximations in general. Usually it is difficult to find a completely satisfactory organi- 
zational principle for a nonperturbative expansion. One that is extensively used is a loop 
expansion. The zero-loop or "tree-approximation" corresponds to the mean-field approxi- 
mation. This approximation is used extensively in field theory and many-body physics. In 
our analysis the first term in Fig. la is identified as the "tree" or mean-field approximation. 
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FIG. 2: The dielectric constant k(/c 2 ) = [1 — ^ + k2 _ m z ] is shown. The singularity at k 2 = 
reflects the operation of the Schwinger mechanism [25]. (Here we see how the infrared singularities 
of the theory lead to nonpropagation of gluons in the QCD vacuum, a result which has long been 
conjectured to be true.) 

That approximation is sufficient to generate the gluon mass via the Schwinger mechanism. 
The second term in Fig. la may be thought of as a (condensate) one-loop correction to the 
mean-field approximation. It is, of course, interesting that we find nonpropagation of gluons 
already in the one-loop approximation. Increasing the number of condensate loops increases 
the number of factors of g 2 0o which appear in the numerator of the terms which make up 
n{k 2 ). That is, at the tree level, we obtain the m 2 G /k 2 term and at the one-loop level, we 
find the term —4r) 2 /(k 2 — m 2 G ). 

We now wish to obtain the contribution to n^(/c) in the Landau gauge of the form 

(5.7) 
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(g^ - k»k u /k 2 ) 



displayed above. Various elements of our analysis are depicted in Figs. 3-5 which are taken 
fromRef. [27]. 

Combining the above result with Eq. (4.3), the mean-field plus the one (condensate) loop 
result for the polarization tensor is 

-Arfk 2 ~ 



rn 



G k 2 -m 2 G 



(51 



We had 



U(k 2 



n(fc 2 ) 

k 2 -r, 



G J 



% / d A xe lk - x < vac | T[j£ a (x)J£ b (0)} | vac> . (5.9) 



From our definition of J®{x), we find 

J a M = gr hc [A»{x)d v Al{x) - 2A»{x)d»Al{x) - Al{x)d»A%x)\ 

+g 2 f abc f a ' b ' c [A b 1 (x)A^x)A a J(x)\. 



(5.10) 



We insert 



(5.11) 



into the last expression to obtain 



= gfa bc [K(x)d»A^x) - 2A^x)d,Al(x) - A^d^x)} (5.12) 
+gf abc [AZ(x)d„A c fl (x) - 2A%(x)d,Al(x) - A°(x)PAfa)] 
+g 2 f abc f a ' b ' c [A b l (x) + Al(x)][A£(x) + A£(x)][Ai(x) + A?(x)}. 

Since we are here working to order (g 2 (j)l), we will drop the last term of Eq. (5.12) at 
this point. (However, we note that it is responsible for the term proportional to m 2 G in Eq. 
(5.8).) Thus, we may use the approximation 



J a u (x) = 9 r bc [A^Al(x) - 2A^At(x) + Kd^Al(x)} 



(5.13) 



for the calculation to be made here. (Note that, in the last term, we have interchanged b 
and c and changed the sign of that term.) We maintain the constraint 



< i) m | d^(x) | i) n >= o, 



(5.14) 
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and implement that constraint by using only the conserved current, Jt^( x ), in our calcula- 
tion. We can define the projection operator 



(5.15) 



which in momentum space has the form 



(5.16) 



We have 



(5.17) 



Thus, using Eq. (5.13), we find 



J a T ,Ax) ~ gf abc [KP^A%x) - 2k p P^d p A p :{x) 



(5.18) 



The first term in Eq. (5.18) is equal to zero, the last term in Eq. (5.18) can be dropped 
because of the constraint 



< i>n I dpA»(x) I i) m >= 0. 



(5.19) 



To the order considered, we have 



% < vac | T[J^ a (x)J^ a ,(x')] | vac > 

= ig 2 f abc f a ' b ' c < vac | T{-2A b p P^d p A c p {x)} x {-2A^P£ V d' p ' A^x')} \ vac > (5.20) 
= i 9 2 f abc f a ' b ' c < vac | A b p I vac > (5.21) 
x{4P7P^ V } < vac | T{d p A c p {x)d' p ' A c p ,{x')} \ vac > 

{AP^P^'}d p d' p ' [iV^,(x,x% (5.22) 



= ig 2 f abc f abc 
We write 



,2 <W fl'p'/ 

" 0O 8 4 



2^ (a;, a/) = S^V^x,^), 



(5.23) 



and find 



i < mc | T[J^ a {x)J^ a ,{x')] | t;ac >= -£ 2 ^cW [AP^P^' d p d' p V^(x,x% (5.24) 
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* v [gluon condensate] 

a b 

\jummimmz/ [ghost propagator] 



[3 gluon proper vertex] 



[4 gluon proper vertex] 



[ghost vertex] 



FIG. 3: Diagrammatic representation of various elements of our model. Note that the dashed line 
contributes for = 0. (See Ref. [27] for further details.) 

Now introduce 



P 1 



+ 



k 2 - m 2 G - U A (k 2 ) (k 2 + ie) 



and note that d p d' p -> [i^] [-ifcj = k 2 . 

Further, P^P^' P T , = P^ v \ and therefore, 



(5.25) 



(WGr'-^'A 2 ) 



n A (fc 2 ) 



(g uu ' - k u k u '/k 2 ) 



-4 V 2 k 2 



k 2 -m 2 G - U{k 2 ) 



<W, (5.26) 



where r/ 2 = (3/32)# 2 
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FIG. 4: Diagrammatic representation of the equations determining the vacuum polarization tensor. 
(SeeRef. [271.) 

n(k 2 ) = mf + TT A (k 2 ) 




FIG. 5: Diagrammatic representation of various elements of the polarization tensor used in this 
work. The first diagram on the right-hand side is responsible for the gluon mass term. [The 
vertex functions in the second two terms are to be expressed in terms of the full gluon propagator, 
%V*(k).\ 

Thus 



Il A (k 2 



[1 



n A (fc2) i 

fc 2 — rri z G J 



21,2 



-Arfk 



_k 2 -m 2 G - U A (k 2 ) 



(5.27) 
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which has the solution, 



n,(^) = (5.28) 

Recall that 

U{k 2 ) = m 2 G + U A {k 2 ), (5.29) 

which then yields Eq. (5.7). 

One may ask how our result is related to the result of a diagrammatic analysis. It may 
be seen that the result given here is obtained if one calculates in the Landau gauge and 
adds a ghost condensate to maintain the transverse structure for U^k). Indeed, Lavelle 
and Schaden [26] have used a ghost condensate to enforce the transverse nature of the 
nonperturbative part of U /xl/ (k). (Recall Eq. (2.26).) Their calculation is made in the deep- 
Euclidean region {k 2 — > — oo). Therefore, we can compare our result with theirs, in the case 
a condensate <: A 2 :> is present, by taking k 2 — > — oo in our result. From Eqs. (5.28) and 
(5.29) we have 

4 

Ii(k 2 ) — > m 2 G - -m 2 G k 2 -> -oo, (5.30) 



3 

m,' 

3 

~32 

which agrees with the result of Ref. [26], when that result is evaluated in the Landau gauge. 
(It is interesting to see how the sign of U(k 2 ) changes as one passes from k 2 = to k 2 = — oo.) 



—ml, (5.31) 
VV)$> (5.32) 



VI. QCD LATTICE CALCULATIONS AND PHENOMENOLOGICAL FORMS 
FOR THE EUCLIDEAN-SPACE GLUON PROPAGATOR 

The form we obtained for the propagator was 

D^{k)=U»- — \D{k). (6.1) 



We now write 



3 k lL —m /S 
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Here Z\ is a normalization parameter which we put equal to 3.82 so that we may obtain a 
continuous representation as we pass from Minkowski to Euclidean space. In Fig. 6 we show 
D{k) with m 2 = 0.25 GeV 2 . (We remark that D(k) = when k 2 = m 2 , D(k) = -Z l /m 2 
at k 2 = 0, and D(k) — > Z x jk 2 for large k 2 .) If we choce Z x = 15.28m 2 = 3.82 our result 
for the propagator will be continuous at k 2 = when we consider both the Euclidean-space 
and Minkowski-space propagators. 

Results for the gluon propagator obtained in a lattice simulation of QCD are given in Ref. 
[28]. In that work the authors also record several phenomenological forms. We reproduce 
these forms in the Appendix for ease of reference. Of these various forms we will make use 
of model A of Ref. [28] which has the form 

AM 2a 



D {k ) 



1 



(k 2 + M 2 Y+ a k 2 + M 2 



L(k 2 ,M) 



with 



L(k 2 ,M) = 



\n(k 2 + M 2 )(k- 2 + AT 



(6.3) 



(6.4) 



and do = 13/22. The parameters used in Ref. [28] to provide a very good fit to the QCD 
lattice data are 



Z = 2Mti, 



(6.5) 



A = 9.841& 



(6.6) 



M = 0.54±1, 



(6.7) 



and 



a 



2.171*, 



(6.8) 



Note that M in GeV units is 1.018 GeV. Rather than work with the lattice data we will 
use Eqs. (6.3)-(6.8) when we compare our results with the lattice data. In Fig. 7 we show 
k 2 D L (k) of Eq. (6.3) and in Fig. 8 we show D L (k). These functions are represented by the 
solid lines in Figs. 7 and 8. Note that Eq. (6.2) may be written in Euclidean space as 

Zi 



D E (k) 



K E + m 3 4+m 2 



(6.9) 
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This form is useful for k\ < 1 GeV 2 and we therefore consider various phenomenological 
forms which may be used to extend Eq. (6.9) so that we may attempt to fit the lattice result 
over a broader momentum range. To that end, we make use of Ref. [29]. The authors of 
that work define the Landau gauge gluon propagator as 

< Al(k)K(k') >= V5(k + k')5 ab (V - ^) (6.10) 

with 

(2 \ 2k 

A 2 k + k2 ) W* 2 )) -7 , (6.11) 

and 7 = —13/22. (We do not ascribe any particular significance to Eq. (6.11). We use 
Eq. (6.11) as a phenomenological form which could be replaced by a form which provides 
a better fit to the data within the context of our model at some future time. We believe 
Eq. (6.11) is useful, since it is a simple matter to remove the first term of that equation and 
introduce a propagator that has the small k 2 behavior of our model.) 

The authors of Ref. [29] introduce two choices for a{k 2 ) of Eq. (6.11). We use their form 
for a 2 (k 2 ): 

Mk 2 ) = — v ^ - xa21 . (6.12) 



In 



e + cii ( 



l QCD , 

In their analysis they put re = 0.5314, A QCD = 354 MeV, a(0) = 2.74, a x = 0.0065 and 
a 2 = 2.40. (Here, we have not recorded the uncertainties in these values which are given in 
Table 2 of Ref. [29].) As we proceed, we will change these values somewhat. As a first step 
we remove the first factor in Eq. (6.11) and write 

Z(k 2 ) = Z 2 (a 2 (k 2 ))^. (6.13) 

We now use ai = 0.0080 and a 2 = 2.10 rather than the values given above. In Fig. 9 we 
show (a 2 (^)) 13 ^ 22 as a function of k, using our modified values of a\ and a 2 . 
We now define 

DEikE) = - 2 2 ))-; m2 . ( a,4) 

K E + m 3 kl+mi 

The function —k 2 D E (k E ) is shown in Fig. 7 as a dotted line. In this calculation we have 
put Z 2 = 2.11. We find a good representation of the lattice result for k E < 2 GeV, 
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1 1 . 1 . 1 . 1 1 1 . 1 

1 2 3 4 5 6 

k 2 (GeV 2 ) 



FIG. 6: The function D(k 2 ) of Eq. (6.2) is shown in Minkowski space. The value for large k 2 is 
given by Z 1 /k 2 with Z 1 = 3.87. Here m = 0.50 GeV. 

In Fig. 8 we compare D E (k E ) with the result of the lattice calculation which is represented 
by the solid line. In Fig. 10 we combine our results in Minkowski and Euclidean space and 
show the values of k 2 D(k 2 ) for both positive and negative k 2 values. For positive k 2 we use 
D(k) of Eq. (6.2) and for negative values of k 2 we use D E (k E ) of Eq. (6.14). Equality of 
these functions at k 2 = implies Z\ = Z2(a(0)) 13 ^ 22 , or Z\ = I.8IZ2. (In our work we have 
used Z x = 3.82 and Z 2 = 2.11. See Eqs. (6.2) and (6.14).) In Fig. 11 we show D(k 2 ) rather 
than k 2 D(k 2 ) ) which was shown in Fig. 10. 

VII. DISCUSSION 

In this work we have developed nonperturbative approximations for the description of the 
gluon condensate and have calculated the form of the gluon propagator. The approximation 
used may be thought of as a condensate-loop expansion. Since the condensate is assumed 
to be in the zero-momentum mode, the loop expansion does not require loop integrals, but 
leads to algebraic relations. Our results are obtained in the Landau gauge. (Note that 
ghosts are introduced to maintain the transverse character of H^k) in Ref. [26].) We are 
able to make some contact with lattice calculations of the gluon propagator, which are made 
in the Landau gauge. We find that our value for the dynamical gluon mass, mg ~ 600 
MeV, is in accord with the results of recent lattice calculations. We have also seen that 
our results agree with those of Lavelle and Schaden, if one evaluates our propagators in the 
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k (GeV) 

FIG. 7: The function —k E DE(k) is shown. The solid line represents the QCD lattice data, while 
the dotted line represents —k 2 E DE{k) in the case that DE{k) is given in Eq. (6.14). 




12 3 4 

k (GeV) 



FIG. 8: The function —DE{k) is shown. The solid line represents the QCD lattice data, while the 
dotted line represents —DE{k) of Eq. (6.14). [See Fig. 7.] 

deep- Euclidean region (k 2 — > — oo) [26]. 

In our work, confinement of quarks and gluons and chiral symmetry breaking are related 
to a single condensate order parameter (g 2 0g). This result is consistent with the fact that 
in lattice simulations of QCD, deconfinement and chiral symmetry restoration take place at 
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0.5 



u.u-| 1 1 1 1 1 1 , 1 1 

12 3 4 

k (GeV) 

FIG. 9: The function (a 2 (k)) 13 / 22 is shown. [See Eq. (6.12).] Note that (a 2 (0)) 13/22 = 1.81. 



4 "I 




k 2 (GeV 2 ) 

FIG. 10: For k 2 > the solid line represents k 2 D(k 2 ) with D(k 2 ) given by Eq. (6.2). Here, 
Z 1 = 3.82. For k 2 < we show k 2 D E (k 2 ), where D E (k 2 E ) is given by Eq. (6.14) with Z 2 = 2.11. 

the same temperature. In this connection, we note that there is no threshold value of (g 2 (fil). 
For any finite value of this parameter, we find chiral symmetry breaking and nonpropagation 
of quarks [30] and gluons. 

In this work we have provided a representation of the gluon propagator in both Euclidean 
and Minkowski space. The Minkowski-space propagator has only complex poles and that 
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k 2 (GeV 2 ) 

FIG. 11: Same as Fig. 10 except that D(k 2 ) is shown. 

implies that the gluon is a nonpropagating mode in the QCD vacuum. Our analysis takes 
into account the important condensate < A^A^ > which is responsible for mass generation 
for the gluon. Our work has some relation to that of Cornwall [31] who obtained a gluon 
mass of 500 ± 200 MeV in his analysis. Cornwall also suggested that "quark confinement 
arises from a vertex condensate supported by a mass gap." 

In recent work, Gracey obtained a pole mass of the gluon of 2.13Ajjg in a two-loop 
renormalization scheme [32]. If we put = 250 MeV, the mass obtained at two- loop 
order in Ref. [32] is 532 MeV, which is close to the value of 500 MeV used in the present 
work. (We remark that in Ref. [22] we obtained a gluon mass of 530 MeV, if we made use 
of Eq. (3.18) of that reference, which includes the effect of including various exchange terms 
in our analysis of the relevant matrix elements.) 



APPENDIX A 



For ease of reference we record various semi-phenomenological forms which are meant to 
represent the Euclidean-space gluon propagator. 
Gribov [33]: 

^(^) = F f^ L( " 2 ' M) - (M) 
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Stingl [34]: 



Marenzoni et al. [35]: 



Cornwall I [31]: 



D L (k 2 ) 



Zk 2 

k 4 + 2A 2 k 2 + M 4 



L(k 2 ,M). 



D L (k 2 " 



Z 



(k 2 ) 1+a + M 2 ' 



D L {k 2 ) = Z 



[k + M (k )] In 



where 



M(k 2 ) = M 



In 



k 2 + AM 2 (k 2 ) 

A2 



-6/11 



A 2 



Cornwall II [36]: 



Cornwall III [36]: 



D L {k 2 ) = Z 



ri2 /fc 2 + 4M 2 \ 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



D L {k 2 ) 



Model A [28]: 



D L {k 2 ) = Z 



k 2 + Ak 2 ln(^) +M 2 ' 



AM 2a 1 , , ' 

+ — — L(k 2 ,M) 



_(k 2 + M 2 y+ a k 2 + M 2 
The parameters for model A are given in Eqs. (6.5)-(6.7). 
Model B [28]: 

D L (k 2 ) = Z 



AM2a 1 r/,2 „^ 



(k 2 ) 1+a + (M 2 ) 1+a k 2 + M 2 ' 



Model C [28]: 



D L {k 2 ) = Z 



A /u2/ir2\a 



M 2 



k 2 + M 2 



L{k 2 ,M) 



(A7) 



(A8) 



(A9) 



(A10) 
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